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ABSTRACT 



Context. Secondary stars in dwarf novae are strongly irradiated during outbursts. It has been argued that this could result in an enhancement of 
^ the mass transfer rate even though the Li region is shadowed from the primary irradiation by the accretion disc. Previous investigations of the 
» I possibility of a circulation flow transporting heat from hot regions to L[ gave opposite answers. 

■ Aims. We investigate numerically the surface flow of irradiated secondaries. We consider the full time dependent problem and we account for 
I the two-dimensional nature of the flow. 
' ' ■ Methods. We use a simple model for the irradiation and the geometry of the secondary star: the irradiation temperature is treated as a free 
parameter and the secondary is replaced by a spherical star with a space-dependent Coriolis force that mimics the effect of the Roche geometry. 
The Euler equations are solved in spherical coordinates with the TVD-MacCormack scheme. 

Results. We show that the Coriolis force leads to the formation of a circulation flow from high latitude region to the close vicinity of the 
Li point. However no heat can be efficiently transported to the L| region due to the rapid radiative cooling of the hot material as it enters the 
equatorial belt shadowed from irradiation. Under the assumption of hydrostatic equilibrium, the Coriolis force could lead to a moderate increase 
of the mass transfer by pushing the gas in the vertical direction in the vicinity of L] , but only during the initial phases of the outburst (about 15 
- 20 orbital periods). It remains however possible that this assumption breaks up due to the strong surface velocity of the flow transiting by Li, 
of the order of the sound speed. In this case however, a three-dimensional approach would then be needed to determine the mass flux leaving 
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Conclusions. We therefore conclude that the Coriolis force does not prevent a flow from the heated regions of the secondary towards the Li 
• • region, at least during the initial phases of an outburst, but the resulting increase of the mass transfer rate is moderate, and it is unlikely to be 
1^ able to account for the duration of lone outbursts. 
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1. Introduction 

Dwarf novae (DN) are cataclysmic variables that undergo out- 
bursts, i.e. a sudden increase of their luminosity by a few mag- 
nitudes lasting for a few days (see e.g. Warner 1995). It is 
now widely believed that these outbursts are due to a ther- 
mal/viscous instability of the disc triggered when hydrogen be- 
comes partially ionized. In the standard model of DN (often 
referred to as the DIM, see Cannizzo 1993; Lasota 2001 for 
reviews), the unstable accretion disc performs a limit cycle be- 
tween an outburst phase of high accretion rate onto the compact 
object and a quiescence phase of low accretion rate where the 
disc re-builds. The role of the irradiation of the secondary star 
during an outburst is still a matter of debate. 

During an outburst the secondary is likely to be signifi- 
cantly affected by the strong increase of the accretion lumi- 
nosity emitted by the primary, as the irradiation flux may ex- 
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ceed by a large amount the intrinsic stellar flux of the sec- 
ondary (see Smak 2004a). It has been often argued, but not 
yet demonstrated, that this could result in a mass transfer en- 
hancement that could, if significant, be an essential ingredient 
of the DIM. For example, it has been suggested that irradiation 
induced mass transfer enhancement is the origin of outburst bi- 
modality (Smak 1999) or that it could play an important role 
in the outburst/superoutburst phenomenon of the SU Uma stars 
(see Hameury 2000; Hameury et al. 2000; Smak 2000). 

However, the Li point is shadowed by the accretion disk 
(Sarna 1990) and thus is not under the direct influence of ir- 
radiation. The very existence of a circulation flow transporting 
heated matter from the irradiated region toward the Lagrange 
point is not clear and one does not know if this results in a 
substantial increase of the mass transfer rate. 

There are some observational clues (see Vogt 1983; Smak 
1995; Smak 2004a) that, in some systems, the mass trans- 
fer rate is increased during an outburst or a superoutburst. 
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Fig. 1. Profile of the Coriolis parameter / along different merid- 
ians of the Roche lobe: = 0° (Li meridian), (p - 5° , 4> - 15° 
and = 45° (from top to bottom). For comparison the profile 
corresponding to spherical geometry is also shown in dashed 
line. 



These observational arguments have however been questioned 
by Osaki & Meyer (2003), who furthermore argue on theo- 
retical grounds against the possibility of such a mass transfer 
enhancement. First, we do not know how the accretion flux 
interacts with the secondary envelope and it is not known if 
irradiation can raise significantly the effective temperature on 
the secondary. Moreover, even in the case where high tempera- 
ture gradients and therefore high pressure gradients are raised, 
Osaki & Meyer (2003) argue that the strong Coriolis force pre- 
vents the formation of a circulation flow transporting heat to- 
ward Li. 

Their main argument is that the Coriolis force remains sig- 
nificant near Li, due to the particular shape of the Roche lobe 
(contrary to the spherical case where the Coriolis force van- 
ishes at the equator); its magnitude remains large enough to 
deflect any meridional flow to what is called a geostrophicflow, 
a steady state where velocity field lines are parallel to isobars 
(see Pedlosky 1982). These conclusions were however chal- 
lenged by computations by Smak (2004a), showing streamlines 
converging to the Li region. In subsequent papers, Osaki & 
Meyer (2004) and Smak (2004b,c) developed their argumenta- 
tion without reaching an agreement. 

One must note that the geometric argument from Osaki & 
Meyer is a one-dimensional argument, valid on the Li meridian 
only. As the L\ point is a singular point in the Roche geome- 
try, the strength of the Coriolis force changes very rapidly as 
one moves away from Li on the equator; this leaves open the 
possibility that the heated material does not follow the shortest 
way along the meridian, but that it follows an "easier" way to 
reach the L\ region. Smak (2004a) computations were done in 
a two-dimensional steady state and thus account for this effect. 

The main flaw of both approaches is that the time scale 
to reach a steady state could be long, possibly longer than the 
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Fig. 2. Profile of the Coriolis parameter / on the Roche lobe at 
colatitudes 6 = 89°, 6 = 85°, 6 = 75° and 6 = 60° (from bottom 
to top). 

duration of an outburst, or that such a steady state could never 
be reached. 

We investigate here the dynamics of the irradiated enve- 
lope by using two-dimensional hydrodynamic simulations. We 
adopt a very simple model for the irradiation and geometry of 
the secondary, but, contrary to Smak (2004a), we solve the full 
time-dependent problem. Our numerical simulations show that 
the flow can indeed reach the vicinity of the Li point, thanks to 
(and not despite) the Coriolis force. However, as hot material 
crosses the boundary between regions directly illuminated by 
the primary star and regions shadowed by the accretion disc, 
it cools down so rapidly that any significant heating of the L\ 
point can not be obtained. If vertical hydrostatic prevails, no 
mass transfer enhancement can occur. However, velocity mag- 
nitudes of the order of the sound speed are obtained in the 
vicinity of the Li point and strong vertical velocity gradients 
are thus expected, which could lead to turbulence and/or the 
breakdown of hydrostatic equilibrium assumption. In this case, 
one would need a full 3D analysis of the problem, far beyond 
the scope of this paper. 

The paper is organized as follow. We present the model in 
Sect.|2]and the numerical code in Sect. [3] Our results are given 
in Sect.|4]and we summarize our conclusions in Sect.|5] 

2. The model 

We consider three models with characteristics of the three DN 
subtypes. The main physical parameters used are given in Table 

m 

2.1. Geometrical considerations on the dynamics 

A complete investigation of the problem would require 3D ra- 
diative hydrodynamic simulations of the stellar envelope in the 
full Roche potential. In order to be tractable, the problem has 
to be simplified. As the scale height of the atmosphere is small 
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Table 1. Physical parameters for the three cases investigated in this paper PqA^ fl> q have their usual meanings. R2 is the mean- 
radius of the secondary. The secondary surface temperature T^, is taken from Smak (2004a). Ti„f) is the maximum irradiation 
temperature reached during an outburst, Y.^^ is the column density of the irradiated layer when T - Ti„. The last column is the 
height of the shadow boundary above the equator when measured on the main meridian. 
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Fig. 3. Map M — > M' from the Roche equipotential to the 
sphere of radius R2 (the mass ratio is ^ = 0.1). The figure shows 
a cut along the Li meridian. 

compared to the stellar radius, <k 1 (even in the vicin- 
ity of the Li point), the vertical dynamical time scale is small 
compared to all other time scales and the vertical hydrostatic 
equilibrium should be reached rapidly. A 2D approach should 
be sufficient to investigate the important features of the surface 
flow (see however the discussion in Sect. l4.3l and in the conclu- 
sion). 

We therefore solve the full set of the Euler equations on the 
surface of the secondary star, whose surface is a Roche equipo- 
tential. As these systems have short orbital periods, the Coriolis 
force is an essential dynamical ingredient. In 2D the Coriolis 
force enters the equations only via its component parallel to 
the surface. Its strength is characterized by the Coriolis param- 
eter f = 2£l.n where n is the local normal to the surface. Fig. 
[U shows the profile of / along dififerent meridians of the Roche 
lobe. On the Li meridian, it can be seen that the Coriolis param- 
eter remains large as stated in Osaki & Meyer (2003, 2004), but 
the 6 dependance of / becomes more and more similar to the 
spherical case as one moves away from Li. Fig. |2]show the 
strong dependence of / on longitude, showing that a ID calcu- 
lation is insufficient. 

The ratio of the inertial force to the Coriolis force is the 
Rossby number: 

V 

9^0=—, (1) 
/L' 



where V and L are characteristic values of the velocity and 
length. In our case the relevant scales are L ~ R2, the mean- 
radius of the secondary, V ~ Cj, the sound speed and f ~ Q., 
yielding Vio ~ 0.04 for all models (see Table[TJ. As expected, 
the Coriolis force dominates the dynamics on the secondary. 

One should in principle use Roche coordinates (see Kopal 
1969) to describe the secondary surface. For simplicity, we 
consider a spherical star of radius R2 (the mean radius of the 
Roche lobe) with a space-dependent Coriolis parameter that 
mimics the Roche geometry. This greatly simplifies the prob- 
lem, as spherical coordinates are much simpler to use than 
Roche coordinates. The mapping between the Roche equipo- 
tential and the sphere is illustrated in Fig.[3j we use a Cartesian 
system of coordinates where the center of mass of the sec- 
ondary is at the origin, the center of mass of the primary is 
at (a, 0, 0) with a the orbital separation and the z axis is per- 
pendicular to the orbital plane. The spherical approximation is 
poor in the Li region, where departure from spherical geom- 
etry is the most important. This has two consequences. First 
the geometrical terms in the Euler momentum equation are not 
well taken into account. However, in this region, the dynamic 
is dominated by the Coriolis force rather than by geometrical 
effect (note that the Rossby number is also the ratio of geomet- 
ric terms to the Coriolis one in the Euler equations). Second 
we underestimate distances: for example the distance from the 
shadow boundary to Li is underestimated by ~ 20 - 30% in 
our model (see Fig. O. This is not a serious difficulty, since it 
is easy to extrapolate our results to a real Roche geometry. 

The Roche geometry is singular at Li, but we assume that 
the stellar surface is at a minimum distance Ar - 0.02/?2 below 
Li, comparable to the vertical scale height of the atmosphere. 
With this assumption, f{6, (p) is no longer singular, as was the 
case in Fig. [T] but instead drops sharply for 9 > 85° and van- 
ishes exactly for 6 - 90° along the L[ meridian (0 = 0°). The 
resulting / profile does not depend much on q (the binary mass 
ratio) and we take here q = 0.1. 

2.2. Irradiation 

It is still debated if the accretion flux can significantly heat the 
secondary atmosphere; we leave this problem for a future in- 
vestigation. Here we model the irradiation of the secondary in 
a crude, but simple way: we consider that the accretion flux is 
completely absorbed and re-radiated as a black body. For con- 
venience, we define Tiir as: 

Tin- = (Fi„it)/cry'' (2) 
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Fig. 4. Time profiles of normal outburst and of superoutburst. 
The superoutburst profile has the same rise and fall but has a 
long plateau at maximum. 

where Fi„ is the irradiation flux and cr is the Stefan-Boltzmann 
constant. Fi„ (and consequently Tiu) has a time dependence 
that account for the temporal profile of the outburst. We de- 
notes by rii-r.o the maximum value of Tin reached during an 
outburst. For simplicity we take Tin.o = 10"^ K for all models 
(see Smak 2004a). We use two outburst profiles, a "short" one 
lasting 60 orbital periods to model "standard" outbursts (note 
that we are not interested in complex behavior such as the out- 
burst bimodality) and a "long" one lasting five times longer in 
order to model superoutburst of the SU Uma subtype (model 
1), see Fig.m Because of the crudeness of the model the exact 
shape and duration of our outburst are not of prime importance. 

Irradiation produces a radiative isothermal layer whose col- 
umn density depends on the irradiation temperature, on top of 
the convective envelope. Using simple vertical structure mod- 
els for the secondary, we have computed the expected column 
density of the isothermal layer when 7^- - 10"* K. The re- 
sults are given in Table[T] The initial column density is small (a 
few tens g.cm"^), because the irradiation flux is negligible, and 
increases with time during an outburst. For the sake of simplic- 
ity, we do not consider the very initial stages of the outburst, 
and assume S - 2^(7" - 10"^ K) as an initial condition. Since 
the precise value of 2 affects only the thermal time, we want to 
obtain the correct order of magnitude of the cooling time scale 
for the gas heated at Tin- - 10'* K. 

The irradiation flux at a given point on the secondary sur- 
face enters as a heating term in the energy equation: 



(3) 



where the intrinsic stellar flux is also included. is a space 
dependent numerical factor accounting for oblique incidence 
and shadowing by the accretion disc on the secondary: 



ip - (n.u) X 5 (el's) 



(4) 



The first effect is accounted for by the n.u term, where n is 
the local normal to the surface and « is a unit vector pointing 



Fig. 5. Profile of if/ along different meridians on the secondary: 
= 0° (Li meridian), = 5°, = 15° and 4> = 45° (from top 
to bottom). The shadow boundary is located at = 80°. Note 
that polar regions are not irradiated. 

towards the primary, assumed to be the source of irradiation. 
The second effect is accounted for by the S (o-j) term, where a, 
is the angle between the orbital plane and the line joining the 
primary center of mass and the running point on the secondary: 



tan Os 



^J(a- xY 



(5) 



-I- 



where {x, y, z) are the Cartesian coordinates of the running point 
(in the same coordinates system as defined previously). We 
then assume that the disc has a fixed opening angle ad (see 
below) and we compute S {as) by: 



5(a,) = ^(l-^tanh^i-^) 
2 Affj 



(6) 



where Aa,/ is an angular size controlling the thickness of the 
transition between the shadowed part and the non-shadowed 
part, taken to be ts.ad - 0.5°. The hyperbolic tangent function 
enables a smooth transition between a, - <k Aa^ and as — 
a 4 » Aaj. 

We use the numerical code from Hameury et al. (1998) to 
compute the vertical structure of the accretion disc during an 
outburst cycle. For each model we determine the maximum 
value of the disc opening angle a^. During an outburst, the 
disc opening is found to be Od ~ 4° for model 1 and ad ~ 6° 
for model 2 and 3. In both cases this translates into a shadow 
boundary located atd ~ 80° when measured on the main merid- 
ian (the larger opening of the disc in model 2 and 3 is compen- 
sated by larger secondaries). The quiescent disc has a much 
lower opening which allows for a much smaller shadowed re- 
gion at the beginning of an outburst. We have checked that 
when the disc enters outbursting state, its opening increases 
very rapidly. Thus it is a good approximation to suppose that 
the disc has its "outbursting thickness" since the very begin- 
ning of the outburst. 
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We use the mapping shown in Fig. [3] to project ^ on the 
sphere of radius R2 in order to obtain ijj - iff{6, 0). Fig.|5]shows 
t// along different meridians. 

2.3. Equations 

We now turn to the mathematical formulation of the equations. 
We use spherical coordinates (6, (p) (6 is the colatitude and (p 
the longitude) to map the sphere of radius ro - Ri- The metric 
properties of these coordinates are given by the metric tensor: 



igu) 



\ rl sin^ I 



(7) 



g'-i denotes the components of the inverse of matrix ^ and g 
denotes its determinant. 

We use the following variables: S is the gas surface density, 
E, H - E + P, and P are the vertically integrated values of 
the energy density, enthalpy density, and pressure, vg, are 
the velocity components. For convenience, we introduce the 
following notations: 

- "Barred" quantities denotes a multiplication by = 

sin 6*. These quantities account for the metric effects; for 
example I^dS - E yfgd9d(p - lA0d(p. 

- "Tilted" velocity components are angular velocity compo- 
nents: Vg - vg/ro = 6 and - v^l{r{) sin0) - (p. 

With these notations, the Euler equations write: 



Pi Pi Pi 

-tvg + — G" + —Gi2^ /tv^+ COS sin 0G22 

dt ^ de dd) ■' two 



d - d - d - 
—E + —Hvg + —Hl't 
dt d0 d(f> ^ 



(8) 
(9) 
(10) 
(11) 



where G'-' - Sv'v-' + Pg'^ is the (symmetric) momentum flux 
tensor. A radiative cooling term is included in the r.h.s. of the 
energy equation in addition to the heating term (O. 

The equations are written in conservative form and differ 
only slightly from the Cartesian case, with the inclusion of 
geometric source terms in the momentum equations (in ad- 
dition to the Coriolis term) and with conservative variables 
(£, £vfl, £v0, E) that are not physical quantities. We also assume 
a perfect gas equation of state: 



P = ZRJ and E = ^ly- + — 
' 2 y- 



1 



(12) 



with an adiabatic index y = 5/3 and Rg - R/fi, R being the 
perfect gas constant. 

Note that we do not include viscosity. As in many astro- 
physical processes, the molecular viscosity is negligible but 
turbulent viscosity could be significant. In a geostrophic state, 
Osaki & Meyer (2003) showed that viscosity is responsible for 
a drift across isobars, but the expected value of turbulent vis- 
cosity yields only a very small drift velocity. 



2.4. Simulation setup 

For each simulation, we start with an isothermal envelope with 
temperature and surface density Sin (see Table [ij. The en- 
velope is initially at rest, i.e. we neglect the quiescence steady 
surface flow feeding mass transfer. The order of magnitude of 
the quiescent flow speed is ~ O.Olcs (see Lubow & Shu 1975), 
small compared to surface velocities found in outburst (see 
Sect. |4]l. Irradiation is then turned on with the temporal pro- 
file simulating an outburst or a superoutburst. Our simulations 
span a time length of 240 PqA- 

3. The numerical code 

The Cartesian-like form of our equations enables us to apply 
directly any scheme developed for Cartesian coordinates. Here 
we use the TVD-MacCormack scheme (see Yee 1987) to solve 
equations (8 - 11) for (L,l,vg,l,v^, E). This scheme is a finite 
difference predictor-corrector scheme followed by a TVD step. 
The "Total Variation Diminishing" step ensures that non phys- 
ical oscillations do not appear (see Hirsch 1990). The two step 
nature of the scheme enables second order accuracy in time as 
well as a second order discretization of the source terms. To ob- 
tain second order in space, the spatial discretization alternates 
between backward and forward differencing during the predic- 
tor and corrector steps. To avoid error accumulation, the order 
of forward/backward differencing is switched at each time-step. 
In steep gradients region, the TVD step adds a diffusive flux lo- 
cally altering the scheme order. The sin cos term in equation 
(6) is discretized in such a way that the physical equilibrium {P 
- est.) is also a solution of the discretized equation. 

Due to the explicit character of the scheme, the numerical 
stability of the scheme is subject to the well known Courant- 
Friedrichs-Levy (CFL) condition; the time step is restricted by 
(see Hirsch 1990): 



Af < Aw = 



1 



VA6»- 



VA0 



(13) 



where c™" - max^Cs/ro + \vg\j and c™"'' - max^Cs/(ro sin0) + 

|v'0|). In our simulations we take Af = /lAfmax with A = 0.95 
(a simulation with A = 0.4 has been checked to give identical 
results). 

We consider a physical domain 0o < < n - 0o and 
-n < (p < 71. We therefore exclude the two polar caps in order 
to avoid difficulties around the poles that are singular points 
in spherical coordinates. This is not a serious problem, as the 
poles are not irradiated (see Fig. |5]). The grid point accumula- 
tion occurring there would severely restrict the CFL condition. 
In the following, we choose 6*0 = 10° and we have checked that 
this particular choice does not influence our results. We use an 
uniform grid with 151 x 200 grid points and we have checked 
the effect of numerical resolution for model 1 on a 76 x 100 
and a 301 x400 grid. As expected the low resolution run shows 
signs of enhanced numerical dissipation and the high resolution 
run shows a more turbulent behavior as compared to our stan- 
dard resolution, but the main features and results presented here 
are not significantly modified by resolution effects (changes are 
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Fig. 6. Snapshots alt - 0.7, 2.5, 5 f orb- Each snapshot shows the temperature field in color gradient (see the color legend), isobars 
and the velocity field. The central panel shows the secondary face on; in the left and right panels, it is rotated by +/- 30 degrees. 
The Li region is marked by a circle whose surface is equal to the stream cross section (see Eq. (17)). For the sake of figure 
readabihty, the magnitude of the velocity field has been multiplied by 3 at f = 0.7 Porb and by 2 at / = 2.5 forb- 



less than 10%). Note that the number of grid points in is odd 
in order to have a grid point at the equator Periodic boundary 
conditions are used at = ±7r and free outflows conditions are 
considered at 6*0 and 7t-0q (we checked that this does not affect 
our results). 

4. Results 

In this section, the velocity, density and temperature at Li refer 
to average quantities computed using a gaussian kernel cen- 



tered on (0 - 90° , (p - 0°) with a FWHM equal to the radius of 
the cross section of the stream leaving the secondary (see Eq. 
dnll below). 

4.1. Model 1: the SU Uma class 

We first discuss normal outbursts. As expected, irradiation in- 
creases the surface temperature of the secondary and, as a re- 
sult, pressure gradients form (see Fig. |6ll. In the presence of 
pressure gradients, the flow tends to reach geostrophism where 
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Fig. 7. Snapshots at f - 18, 30, 45 forb- Each snapshot shows the temperature field in color gradient (see the color legend), isobars 
and the velocity field 
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Fig. 8. Left panel: Evolution of the velocity at Li normalized by the isothermal sound speed (computed with the temperature at 
Li) during an outburst (solid line) and a superoutburst (dashed line). The maximum velocity is reached at f = 18 Porb (see Fig. 
|7]i. Right panel: Evolution of the surface density at the Li point normaUzed by its initial value surface density during an outburst 
(solid line) and a superoutburst (dashed line). 



the Coriolis force balances pressure gradients (i.e. velocity is 
along constant pressure line). This dynamical phenomenon is 
known as "Rossby adjustment" or "geostrophic adjustment" 
(see Pedlosky 2003). Initially, vertical pressure gradients drive 
a vertical acceleration of the fluid, pushing it towards the equa- 
tor However the Coriolis force quickly deflects the flow west- 
wards and this effect is stronger in the vicinity of Li. This 
leads to an oscillatory pattern of the velocity field (see snap- 
shots at f = 0.7 Porb in Fig. |6]l along the shadow boundary of 
the disc. This pattern drives the flow toward Li. In addition, 
vertical pressure gradients appearing at high longitude, where 
the normal irradiation flux drops, are also adjusted and as a re- 
sult, a clockwise circulation flow sets in. At t - 2.5 Porb (see 
Fig. |6]l, the flow has reached geostrophism. The flow crosses 
isobars only in the close vicinity of Li, as the Coriolis force is 
weaker there (see Sect. 12. II ) and the flow speed larger Note that 
the geostrophic state is significantly different from the one ex- 
pected from the initial pressure gradients; this shows the com- 
plexity of the geostrophic adjustment process during which 
both the flow velocity and the pressure gradients vary, the ir- 
radiation flux being the only fixed quantity. 

The overall situation can be depicted as follows: irradiation 
drives a giant anti-cyclonic perturbation on each hemisphere of 
the secondary. Furthermore, the associated clockwise circula- 
tion goes through the Li point. In the same time, an equatorial 

jet traveling westwards forms at longitude <p 60° as well as 

a weaker eastwards one traveling one at ~ 60° (see snapshots 
at t = 2.5 Porb in Fig.|6]l. 

As time goes on, pressure gradients increase with irradia- 
tion and, as a result, the circulation strengthens (see Fig. [7] at 
? = 30 Porb; the outburst maximum). The left panel of Fig. |8] 
shows the magnitude of the velocity at Li . The maximum ve- 
locity is reached at f = 18 Porb where ~ 0.6cs (for symmetry 
reason vg - 0). The corresponding snapshot is shown in Fig.|6] 
The right panel of Fig.[8]shows the evolution of the surface den- 



sity at Li . In the early stage of the outburst, the surface density 
increases by a few percent as matter is initially driven toward 
the equator Later on, it decreases temporarily by a few percents 
due to the strong surface flow through Li and as a consequence 
of mass conservation. 

Snapshots at different times (see Figs. |6] |7]l show that the 
anti-cyclonic perturbation slowly drifts westward. This is a nat- 
ural consequence of its clockwise circulation and this shows 
that no steady state can ever be reached. As a consequence, 
the circulation flow initially transiting by the Li point moves 
away from Li (see left panel of Fig. |8]for f < 50 Porb)- Note 
also that the temperature distribution at the surface of the sec- 
ondary is not much aff'ected by the flow, as would be the case if 
heat transport by advection were important. This is an impor- 
tant point that will be discussed in more details later on. 

The subsequent increase of the velocity at Li in Fig.|8](left 
panel, at f ~ 70 Porb) is due to the transit of the westwards jet. 

At the outburst end, the residual circulation flow slowly de- 
cays (due to numerical dissipation) on a time scale shorter than 
the time interval between outbursts; physical viscosity (e.g. due 
to turbulence) would bring the velocity field to zero before the 
next outburst starts. 

Numerical simulations of a superoutburst shows the same 
features during the initial phase, as expected. Even after t - 
30 Porb; when a normal outburst would have started to decline, 
the surface flow is not very different in the normal and super- 
outburst cases, see Fig. |7] Due to the maintained irradiation, 
the strength of the surface flow does not decay as it does in the 
normal outburst case. However, the anti-cyclonic perturbation 
drifts westward in the same way as it does in the normal out- 
burst case, also suppressing the flow transiting by Li (see left 
panel of Fig. [8]l. The subsequent evolution shows a turbulent 
like behavior that slowly decays when the superoutburst ends. 
The surface density at Li (see right panel of Fig.[8]is increased 
by at most 20% during a superoutburst. 
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Fig. 9. Left panel: Particles pathlines from f = to f = 40 Porb- Right panel: Time evolution of the colatitude of each particles. 
The shadow boundary is at ^? = 80°. 



In order to compare our results with those from Smak 
(2004a), we now follow the trajectories of test particles dis- 
tributed in the atmosphere and initially at rest. Fig.|9](left panel) 
shows a few example of such trajectories. Each trajectory lasts 
for 40 orbital periods. This figure shows that particles initially 
at longitude (p ~ 40° - 50° reach the close vicinity of L]. The 
arrival time in the Li region can be obtained from the right 
panel of Fig. |9] which shows the time evolution of each parti- 
cle colatitude. Fig. [10] shows the temperature of each particle 
vs colatitude. The particles are heated during their trajectory 
towards L\ and they cool down very quickly once they enter 
the shadowed region. As a consequence, the temperature of the 
Li point remains very close to its quiescent value (2500 K), 
contrary to Smak (2004a). The reason for this is that the cool- 
ing time scale is shorter than the transit time towards Li: as 
shown in the right panel of Fig. |9] the time scale for crossing 
the shadow boundary is of the order of few orbital periods. The 
cooling time scale is: 

Tcooi = -ir (14) 
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Fig. 10. Temperature vs colatitude for each particle (To = 2500 
K, the shadow boundary is at = 80°). Particles cool very 
rapidly as they enter the shaded region. 



where c,. is the specific heat capacity. Thus TcooI °^ with 
a very strong temperature dependance. For an initial value Z = 

350 g.Cm-2, Teool = 9 Porb for T = 2500 K but Teool = 0.15 Port 

for T - 10"* K. Consequently as materials heated at T ~ 10"^ K 
enters the shadowed region, their very rapid cooling is unavoid- 
able. A more elaborate computation of a T'* cooling shows that 
after a few cooling time scales the gas temperature is close to 
the ambient value. The preceding argument holds even better 
when one corrects for the Roche geometry (see Sect. |2]l: real 
distances are about 30% larger. The cooling of the gas is so ef- 
ficient that this conclusion does not depend much on the thick- 
ness of the shadowed region (see discussion at the end of Sect. 
I2.2I 1 nor on the exact value of Sin, since it would require un- 
realistically large values of Z for TcooI to be of the order of the 
crossing time scale. 



4.2. Model 2 & 3: U Gem and Z Cam classes 

Numerical simulations with parameters corresponding to U 
Gem (model 2) and Z Cam (model 3) show very similar fea- 
tures to the SU Uma case. In particular the time needed for 
particles to cross the shadow boundary is of the same order as 
previously found, namely a few orbital periods while the cool- 
ing time is unchanged: for the parameters of model 2 and 3 
Teool ~ 0.02 - 0.03 f 01b (due to larger orbital periods), thus 
showing that the above conclusion also applies here. 

4.3. Mass transfer enhancement ? 

We now discuss the implications of our numerical results on a 
possible mass transfer enhancement. The mass transfer rate is 
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Fig. 11. Vertical profile of the density p in the isothermal envelope when vg -Q,vg - 0.5cs, vg - 0.75cs and vq — (from bottom 
to top), left panel corresponds to Z£, - 1.2H and right panel corresponds to zl, - 2.H. The density is normalized by the value of 
the density at Li when vy = 0. 



given by (see Lubow & Shu 1975, Meyer & Meyer-Hofmeister 
1983): 



M= ep(L,)Cs 

where Cs is the isothermal sound speed: 



(15) 



(16) 



and Q is the cross section of the stream (see Meyer & Meyer- 
Hofmeister 1983): 



Q 



27r/Cs\2 

T'q'' 



(17) 



where ^ is a constant which depends only slightly on the mass 
ratio q. 

The factor Qc^ yields a T{L\fl'^ dependence of the mass 
transfer on the temperature at Li. If we further assume that 
vertical hydrostatic equilibrium prevails, p(Li) also depends 
T{Li) and on 2(Li). Our results have shown that T{Li) is al- 
most unaffected by uradiation and heat transport. An increase 
of S(Li) by a few percents leads to an increase of M by the 
same amount and is therefore inefficient to increase the mass 
transfer rate (see right panel of Fig. [8]). However, beneath L\ 
the Coriolis force pushes the gas in the vertical direction and 
this could change the vertical distribution of the mass density. 
Including the Coriolis force in the hydrostatic balance gives: 



az 



(18) 



where the z axis now refers to the local vertical axis with the 
origin taken at the base of the isothermal atmosphere, g is the 
Roche gravity and vy is the surface velocity of the gas. 

Due to the Coriolis force, the location where the effective 
gravity vanishes can be lower than the L\ point. In this case, 
one should replace p{L\) in Eq. [T5]by the density at the point 
where the effective gravity vanishes. 



When V|| — 0, Eq.[T8]can be solved by expanding g to first 
order around zli , the height of Li (see e.g. Lubow and Shu 
1975), giving: 



po(z)=Po(ii)exp((z-Zi,)2///2) 



(19) 



where H - Cs /(Q + 1 /2) is the vertical scale height of the 
isothermal envelope (A is a numerical factor given in Lubow 
and Shu (1975)) 

Since we do not known the vertical profile of vy, we solve 
Eq. [18] with vg = 0.5,0.75, 1.0 Cj. The results are shown in 
Fig. [TT] The mass transfer enhancement is directly given by 
the increase of the density. Results depends on the ratio zli /H, 
which has to be determined by computing the vertical structure 
of the envelope. This is beyond the scope of the present paper 
and we consider in Fig. [TT|two cases: ZlJH - 1.2 wdziJH - 
2. 

It should be noted that, as the velocity reaches a fraction 
of the sound speed, it is possible that the hydrostatic equilib- 
rium assumption breaks in the whole vertical extend of the en- 
velop. A detailed analysis would require 3D simulations, out of 
the scope of this paper However, the main conclusion that the 
cooling time is much shorter than the crossing time of the disc 
shadow will not change. 

Finally, note that a mass transfer enhancement due to the 
strong surface flow, whatever the mechanism at play, would be 
restricted to the rise of the outburst or superoutburst (see left 
panel of Fig.lfor t < 30 Porb)- 

5. Conclusion 

We have numerically investigated the surface flow of the irra- 
diated secondary star in dwarf novae, using a simple model for 
both the geometry and irradiation of the secondary star. 

Our numerical simulations can be viewed as large scale me- 
teorological simulations of irradiated secondary stars: irradia- 
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tion triggers a large anti-cyclonic perturbation (i.e. high pres- 
sure perturbation) with the associated clockwise circulation 
due to the Coriohs force. An important result of our numerical 
simulations is that the Roche geometry naturally leads to a cir- 
culation flow transiting by Li, contrary to naive expectations. 
This provides a way to transport matter/heat in the vicinity of 
Li, where matter leaves the secondary. We have also shown 
that it is necessary to solve the full time-dependent problem 
in order to catch the dynamics of the geostrophic adjustment 
process; the flow does not reach the geostrophic state without 
modifying the background pressure gradients (as is implicitly 
done in Osaki & Meyer 2003; 2004). It is also important to 
note that no steady state can be reached on the time scale of an 
outburst or a superoutburst. 

We modeled the irradiation and radiative cooling of the en- 
velope of the secondary in a very crude way, thus avoiding the 
necessity of solving the vertical structure of the atmosphere. In 
our model, radiative coohng occurs in a very short time scale, 
much shorter than the crossing time scale of the shadowed re- 
gion. There is thus no efficient heat transport in the vicinity 
of Li . It is unhkely that this conclusion will change in a 3D 
approach that would also include a more realistic treatment of 
radiative transfer. As a consequence, if hydrostatic equilibrium 
prevails, the only effect that could lead to a moderate enhance- 
ment of mass transfer is the vertical Coriohs force. It could 
however be possible that hydrostatic equilibrium breaks down 
as the magnitude of the surface flow becomes of the order of the 
sound speed. A 3D approach would then be necessary to deter- 
mine the mass flux leaving the secondary. Note that we have 
not considered here the magnetic field of the secondary. For 
the expected range of intensity, the magnetic field is unlikely to 
dominate the gas dynamics but its influence on the flow could 
be of some importance. An correct implementation of the mag- 
netic field configuration would however be very difficult and is 
far beyond the scope of the paper. 

Apart from the transport mechanism investigated here, any 
other mechanisms that could directly heat the Li point would 
also contribute to an increase of the mass transfer rate. For ex- 
ample the direct heating of the Li point by the rim of the ac- 
cretion disc or by the scattering of the accretion luminosity by 
outflowing matter are currently under investigation to check if 
heating of L\ could be efficient. We leave this for a forthcoming 
paper. 

Even if the circulation flow found here were to lead to a 
significant increase of the mass transfer rate, this effect would 
not last longer that 10 - 20 orbital periods, as the anti-cyclonic 
perturbation moves rapidly westwards, quenching the flow 
through Li on a time scale comparable to the rise time of the 
outburst. We would expect in this case no difl'erences between 
outbursts and superoutbursts. 

Finally, it could be tempting to apply these results to soft X- 
ray transients (see e.g. Chen et al. 1997 for a review). In these 
systems the irradiation flux is larger than in DN by a factor of 
order 10^ and consequently 7"^ ~ 5.10"' K. Such a strong irra- 
diation raises very important pressure gradients and one could 
expect supersonic velocities. A numerical investigation of this 
case would need a more robust numerical code. 
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